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Abstract 

We introduce new 6D standing wave braneworld model generated by gravity coupled to 
a phantom-like scalar field and investigate the problem of pure gravitational localization 
of matter fields. We show that in the case of increasing warp factor spin 0, 1, 1/2 and 2 
fields are localized on the brane. 

PACS numbers: 04.50.-h, ll.25.-w, 11.27.+d 



1 Introduction 

The scenario where our world is associated with a brane embedded in a higher dimensional 
spacetime with non-factorizable geometry has attracted a lot of interest with the aim of solving 
several open questions in modern physics (see [Hll] for reviews). 

The braneworld models assume that all matter fields are localized on the brane, whereas 
gravity can propagate in the extra dimensions. Recently we investigated the standing wave 
braneworld model in 5D spacetime [5H9] and had shown that it provides universal gravitational 
trapping of zero modes of all kinds of matter fields in the case of rapid oscillations of standing 
waves in the bulk. The goal of this article is to generalize the model to 6D spacetime. 

Although there exist a vast literature concerned to the localization of fields in 6D braneworld, 
there is not yet found a universal trapping mechanism for all fields. In the existing 6D models 
with stationary exponentially warped spacetimes spin-0, spin-1 and spin-2 fields are localized 
on the brane with the decreasing warp factor, but spin-1/2 fields can be localized only with the 
increasing warp factor [10l[Tl]. There exist also 6D models with non-exponential warp factors 
providing gravitational localization of all kinds of bulk fields on the brane [T211T9] . however. 
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these models require introduction of unnatural gravitational sources. There have been also 
considered models with time-dependent metrics and fields [SOHSS]- 

Here we introduce the non-stationary 6D braneworld model with gravity coupled to a 
phantom-like scalar field in the bulk, where generated bulk standing waves are bounded by 
the brane at the 2D extra space origin and the static part of the gravitational potential in- 
creasing at the extra space infinity. Then we explicitly show that this model provides universal 
gravitational trapping of zero modes of all kinds of matter fields in the case of rapid oscillations 
of bulk standing waves. It must be mentioned that analogous physical setup was considered in 
recent paper [21], which differs from our model in metric ansatz. 

In the article in Section 2 we introduce the background solution of our model and some 
expressions for using in subsequent sections. Then, in Sections 3, 4, 5 and 6 we demonstrate 
existence of normalizable zero modes of spin-0, -1, -1/2 and -2 particles on the brane. Short 
final conclusions can be found in Section 7. 



2 Background Solution 

In 6D spacetime the Einstein equations with a bulk cosmological constant A and stress-energy 
tensor Tab are 

Rab - -gAsR = J^gAB + k'^TAB, (1) 

where capital Latin indices refer to 6D spacetime and 6D gravitational constant k obeys the 
relation k"^ = SttG = {G and M are the 6D Newton constant and the 6D Planck mass scale, 
respectively). For the two extra spatial dimensions we introduce polar coordinates {r,6), where 
< r < cx) and < ^ < 27r. 

In our physical setup there exists non-self-interacting phantom-like scalar field, which de- 
pends on time, propagates in the bulk and has the following energy-momentum tensor 

Tscal.fleldAB = -C?A</'C?iJ</> + ^gABd'^(f)dc(f) ■ (2) 

Using ([2]), equations ([T]) can be rewritten in the form 

Rab = -^AgAB - k^dA(pdB(p ■ (3) 
We look for the metric in the form 

ds^ = e^- [dt' - e"(*'^) {dx^ + dy' + dz')] - dr^ - i^^e^— 3n(i,r)^^2^ 

where curvature scale a 7^ and i?o > are constants, and function u depends only on time t 
and on extra radial polar coordinate r. Then u{t,r) = this metric ansatz coincides with the 
known solution [11] , that describes the model with a 3-brane at the origin r = which is a 4D 
local string-like topological defect in the 6D spacetime. We also mention that the metric ansatz 
(jlj) differs from that recently proposed in the article [21], our metric symmetrically considers 
X, y and z coordinates. 
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Taking into account symmetry properties of the metric ansatz @ we assume that phantom- 
hke scalar field depends only on time t and extra coordinate r, i.e. (p = (j){t,r). After a 
straightforward calculation, equations reduce to 



^2ar _ l\^2ar j_ l.2d^4> 

odudu _ i.2d±d± 
dt dr at dr ' 

- lOa^e^^'^ - hae^"'-^ - e^^^-f-l = Ae^"^ (5) 



ot^ or or-^ 



i^i 

from which we get 



A = -10a' 



"^Kat) ~ \dt) ^ (P.\ 
r^dudu ^ j^2d±d± 
pt dr dt dr'' 

d u c:„„2ar&u ^2ar d u n 

dt^ e - u. 

First equation in system ([6]) fixes relation between bulk cosmological constant A and the 
curvature scale a in the exponential warp factor of the metric (j4]). We see that the bulk 
cosmological constant A must be negative. 

Taking into account the second and the third equations in ([S]), we set the relation between 
metric function u {t, r) and phantom-like scalar field {t, r) as 

<Pit^r) = ^u{t,r). (7) 

Now we return to the last equation in dH]) and, using separation of variables u {t, r) = 
T (t) f (r), decouple it as follows 

r + cjV = 0, /" + 5a f + u^e-^"' f = 0, (8) 

where where > is some real constant and overdots and primes denote derivatives with 
respect to t and r, respectively. The solution to the first equation in is 

r (t) = ci sin (ut) + C2 cos (ut) , (9) 

where ci and C2 are some real constants. To solve the second equation in (|8]) we perform the 
following change of variables 

z = ^e-^^, f{r) = zlq{z), (10) 



and it gets the form 



with the solution 



dz'^ zdz \ z 



q{z) = ArJs{z) + A2Y,{z), (12) 
3 



where Ai and A2 are some real constants, and Js (2) and Ys (z) are |-order Bessel functions 
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of first and second kind respectively. Taking into account change of variables (fTOl) . for /(r) we 
get 

/ (r) = Cse-i'^Vs (^^"'^'^^ + C4e-i'^"F5 (^e""''^ , (13) 
where C3 and C4 are some real constants. So, in general, u{t,r) has the form 

M(t, r) = (ci sin (wt) + C2 cos (ut)) (cse-l^'J^ f ^e-"" ) + c^e-^'^'Ys ( ^e-^A ] . (14) 
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In what follows we will consider the case of increasing warp factor, i.e. a > and choose 
M(t, r) in the form 

u{t,r) = sin {ut)B{r) , (15) 

with 

B{r) = Boe-^^'Ys f-e""'') , (16) 

2 V (2 / 

where Bo is some real constant. 

The ghost-like field (j){t, r) and the metric oscillations via the metric oscillatory function 
u{t,r) must be unobservable on the brane located at r = in the extra space. According to 
(El), ( fT5|) and ( fT6l) we can fulfil the requirement by imposing the boundary condition on the 
Bessel function 

Fs f-e""'"') =0. (17) 

2 Va / r=0 

Taking into account the oscillatory character of the Bessel function Ys and denoting by Z„ its 
n-th zero, this condition can be written in the form 

- = Zn, (18) 

a 

which quantizes the standing wave oscillation frequency u in terms or curvature scale a. In our 
case of increasing warp factor, i.e. a > 0, and for the some chosen fixed value of n in (ITS!) . in 
the 2D extra space there can exist finite number of circles where the oscillatory function Ys 
vanishes. The radiuses of such "zero circles", which correspond to the nodes of the standing 
waves, satisfy the relation 

-e-"'^ = , (19) 
a 

where i = 0,1,.. .,n — 1. This "zero circles" can be considered as the islands where matter 
particles can be bound in the extra space. Obviously the metric function u(t, r) f ll6p and scalar 
field <^(t, r) ([7]) will be also vanishing at the "zero circles". In what follows in f|T8|) we assume 

n=l, i.e. we choose the first zero of Ys, which is equal Zi ^ 3.96. In this case standing waves 

2 

have only one node at r = 0. 

In the matter field equations the metric oscillatory function f|T5|) . oscillating with the fre- 
quency of the standing waves, enters via exponents e*"*^*''') with b denoting some real constant. 
We assume that the standing wave frequency u is much larger than the frequencies correspond- 
ing to the energies of the particles localized on the string-like brane. Under such conditions in 
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matter field equations one can perform time averaging of oscillating exponents. Denoting by 
Kf,{r) the time average <^e^''(*''')) and using the results of our previous papers [HlElin], we can 
write the following useful expressions for time averages 



H-(t>-(t>-(^) 



= , 



(20) 



where Iq is the modified Bessel function of order zero and the function B[r) is defined by f|T6|) . 
At the origin and infinity function Kh{r) has the following asymptotic forms 



2TT{uj/a) \ ' 

T ( 3^2 bB(i \ I ^ (p-a.r\ 



(21) 



3 Localization of Scalar Fields 

We start with the problem of localization of massless scalar fields defined by the 6D action 
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where the determinant for our ansatz (HI) has the following form 



(22) 



-g = Roe 



5ar 



(23) 



Performing time averaging of the oscillating functions in the corresponding Klein-Gordon equa- 
tion, 

^ =5m(v^/'^5^$)=0, (24) 



we get 
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- K, {dl + dl + d^:) - ^,K,dl 



(25) 



where functions Ki{r) and K^{r) are time averages defined by fl20|) . We look for the solution 
to this equation in the form 



$ {t,x,y,z,r, 



J(Et-pxX-pyy- 



l,m 



(26) 



with (pmif) obeying the following equation 



(r) = 



(27) 
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where = vi + + vl- the brane, where m ~ 0, the parameters E, p^, py and Pz can be 
considered as energy and momentum components along 4D. 

Taking into account that the energy and momentum of 4D massless scalar field obey the 
dispersion relation 

E^=pI+pI+pI, (28) 
for 00 (r) of the s-wave (/ = 0) zero mode wave function we get the equation 



r ^2 



d 



Mr) = 



(29) 



We put ( 129|) into the form of an analogue non-relativistic quantum mechanical problem by 
making the change 

(30) 
(31) 
(32) 



and for <Jo{r) we find 



where the function 
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ao = 



Uoir) = ^a'+p'e-"'^IU-l) 



is the analog of non-relativistic potential. Figure 1 shows behaviour of Uq^t). 



14- 



12- 



10- 
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Figure 1: The effective bulk potential (l32l) . 
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We explore fl3Tl) in two limiting regions, far from and close to the brane. Using f l2T|) . far 
from the brane, i.e. r — )■ +00, the equation flHT]) obtains the form 



r 25 2' 



c^o = , (33) 



with the general solution 



(r) = Cie-^'^''/^ + C2e^"'■/^ (34) 

where Ci and C2 are some constants. Taking into account (l30l) . the extra dimension factor 
</)o(r) of the scalar field gets the following form 

<Po{r)l^+^ = C,e-'^^ + C,. (35) 

To have normalizable zero mode, at the infinity we impose the boundary condition 

</'o(r)|,^+^ = 0. (36) 

So we set C2 = in ( l35l) and get the following result 

0o(r)|,_+^ = Cie-^- . (37) 

Taking into account (12T|) . close to the brane, r — 0, the equation (!3T|) again reduces to (!33|) . 
Resemblance of the equations far from and close to the brane is not surprising since according 
to ( I2TI) the second term of U (r) in ( l32l) vanishes in both regions (these two regions correspond 
to the nodes of the standing waves). Therefore for the extra dimension factor (f)o{r) of the scalar 
field zero mode we again get 

0o(r)|,_o = C3e-^'^^ + C4, (38) 
where C3 and C4 are some constants. At the origin r = we impose boundary condition 

d(j)o (r) 



dr 



= , (39) 

r=0 



and get the following result 

0o(r)|,^o = C3(5a + e-^'^^) . (40) 

So 4>o{f) has a maximum on the brane and falls off at the infinity as e~^°'^ . In the scalar 
field action ( |22l) the determinant ( |23l) and the metric tensor with upper indices give the total 
exponential factor e^"^, which obviously increases for a > 0. But the extra part of wave function 
(13 7p contains the exponentially decreasing factor e~^"^. So, for such extra dimension factor the 
integral over r in the action fl22l) is convergent, i.e. 4D scalar fields are localized on the brane. 
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4 Localization of Vector Fields 

For simplicity we consider only the U{1) vector field, the generalization to the case of non- 
Abelian gauge fields is straightforward. The action of vector field is: 

S, = --! d'x^ g''''g''''FMpFNR , (41) 



where 

Fmp = OmAp - dpAM (42) 

is the 6D vector field tensor, and the determinant g of the metric is defined by (12^ . 
The action (jH]) gives the system of six equations 

9,/ (v^ = , (43) 



which for our metric (jl]) have the following explicit form 



e" {d^F^t + dyFyt + d,F,t) + dr (e^-F,^) + ^^e'^^'+^^deFt 



et 



dt ie"Ft.,) - dy (e'^'-^^Fy,) - (e'^^-^^F,,) - dr (e^— "F^) - ^^e^r+2uQ^p^^ ^ g, 
dt (e^^-^Fty) - {e^^-^^F,y) - (e— ^"F.j,) - 9, (e^— - ^^e'^^+^^deFey = 0, 
dt {e"Ft,) - (e»'-2"F,,) - dy (e^^'^^F,,) - dr (e^^-^F,,) - ^e'^^+^^a.Fe, = 0, 
dt (e^-Ft,) - d., (e^— "F„,) - dy (e^^Fyr) - d, (e^^F.r) - ^e^'^'+'^^deFer = 0, 
dt {e^'+^'^Fte) - d^ (e"^+2"F,,) - dy (e^^+^^Fye) - d, (e^'^+^^F,,) - dr (e^^+^^F.^) = 0. 

We look for the solution to the system fHS]) in the form: 

^(x^) = atix-')Zpmir)e''', Af,{x^) = e"(*'^)afc(x-) E Pm(r)e^'^ 

l,m 1,171 

Ar{x^) = a,(x'^) Eem(r)e^'^ ^(x^) = , 



(44) 



(45) 



Lm 



where index k runs x, y and z, u{t,r) is the oscillatory metric function ( !T5|) . and a^(x'^) denote 
the components of 4D vector potential (Greek letters are used for 4D indices) that obey 4D 
Lorenz-like gauge condition 

7]''^daaisix'') = (46) 

(?7q/3 denote the metric of 4D Minkowski spacetime). In fact, equation fj46|l together with the 
last expression in psjl can be considered as the full set of gauge conditions imposed on the 
vector field Am{x'-^)- 

In the case when the frequencies of bulk standing waves u is much larger than frequencies 
associated with the energies of the particles on the brane, 

w > F , (47) 
we assume the existence of localized flat 4D vector waves, 

{x") ~ £^eKi?i-P.^-P«s/-P.^) ^ (48) 
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where p^, py, p^ are components of energy-momentum along the brane. So putting fjlSj) into 
and performing time averaging we get the following system of equations 



(49) 
(50) 



a^a^ + (Ki - 1) dt + -2 ) = , (51) 



dr (e^-e^) = , (52) 

where in (l50i) index i runs values x, ?/ and ^, and the functions ii'i(r), -ft'3(r) and K4^{r) are 
time averages defined by (!20|) . For the s-wave (/ = 0) zero mode wave function this system 
reduces to 

{dr (e3-9,) - eV {K, - 1)) a^po = 5. (e^'^'^eo) -9<a, , (53) 
{dr (e3-9,) - eV (Ki - 1)) a,po = i^i^r (e^^'-^o) d,ar , (54) 

p2(A'i-l)a,eo = , (55) 
dr (e^-eo) = , (56) 

with p^ = Px+pI +^2- According to (!55l) and (156!) we set ^o('^) = Cie~^°''^ and ar(x'^) = 0, then 
equations (l53l) and (154|) reduce to the following single equation 



So the zero mode solution has the following 6D form 

A{x^)=atixnPoir), A,(x^) = e"(*'^)a,(a:>o(r), A,(x^) = 0, A,(x^) = 0, (58) 
where 4D factors a^(x'^) have form defined by (HHl) . and the function po(?") obeys the equation 
By making the change 

Po(r) = e-3-/2^o(r), (59) 
we rewrite (1571) into the form of non-relativistic quantum mechanical problem 



Vo = 0, (60) 



where the potential 

i7i(r) = ^a2 + pV2-(i^'i-l) (61) 
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differs from the analogous potential for scalar fields fl32l) (see Figure 1) only in constant factor 
in the first term. So, using the same arguments as in the previous section, the solutions to f lBU]) 
in the two limiting regions - far from and close to the brane - are 



vo (r 



Vo (r) 



(62) 
(63) 



where Ci, C2, C3 and C4 are again some constants. Taking into account (!59l) and imposing 
boundary conditions analogues to (136|) and (139|1 . the extra dimension factor po{r) of the vector 
field zero mode will have the following asymptotic forms 



Po [r) 



Cie- 



-3ar 



Po [r) 



Cs{3c 



-3ar 



) 



(64) 



So po{r) has maximum on the brane and falls off at the infinity as e ^"'^ . Now, using 
and P2|) . it's easy to find the time-averaged components of the zero mode tensor 



dr 



k 

i ) 



(65) 



where f^u{x") = 9^a^(x") — 9ya^(x") and small Latin indexes i and /c run values x, y and z. 
Accordingly the 6D Lagrangian for zero mode is 



(0) 



^{F^) (Fi) = -f [po'fUfi: + Po' (K, - 1) ifid'a, + ^d^a,) + 

+ po\Ki - ifd'atd'a, + (^)'a^a'^ 



(66) 



Using ( 12T|) and (l64l) . it's straightforward to show that on the brane, r = 0, the Lagrangian has 
standard 4D form 



r=0 



-fC3^(3a + l)V;/.^. 



and far from the brane, r — 00, it has the following asymptotic form 



(0) 



/~i 'A ~5ar fv fa 
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(67) 



(68) 



So, the integral over extra coordinates r and 6 in the action (HTi) is convergent, what means 
that 4D vector fields are localized on the brane. 



5 Localization of Spin 1/2 Fermionic Fields 

The purpose of this section is to explicitly show that there exists normalizable zero mode of 
the spin 1/2 fermionic field. The starting action is the Dirac action 

Si=f d^x^/^^lT'^DA'^ , (69) 
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from which the equation of motion is 

= r^D^^ = (r^D^ + + v^Dg) ^ . (ro) 

We introduce the vielbein through the conventional definition: 

9ab = VAshih'^ , (71) 

where A,B, refer to 6D local Lorentz (tangent) frame. Using the relation = h^'j^ with 
and 'j^ being the curved gamma matrices and the flat gamma ones, respectively, we have 
relations 

r* = e-""7*, = e-"/^Y, r = Y, = i?o ^e-'^"+=^"/%^ (72) 
where the index i runs values x, y and z. The covariant derivatives in ( l69i) and ( 170|) are 

DA = dA + \ni^lBlc , (73) 

with Q^'-' being spin-connections. The nonvanishing spin-connection components for our back- 
ground metric @ are 

Taking into account these results the Dirac equation can be written as 

+ {^d, + Ydy + Yd,) + e'^'Y ( Y + + ^^^^l^'-^^,) * = . (75) 
We look for solutions of the form 

v^(x^)=^(x^)^«^(r)e^'^ (76) 

l,m 

where ■ip{x'^) satisfies the massless 4D Dirac equation 7^*9^?/; (x*^) = 0. Using f l76|) . for the s-wave 
(/ = 0) zero mode fermionic wave function Dirac equation fl75|) reduces to 



1 {Yd, + Ydy + Ydz) ^A(a;^)ao(r) = . (77) 



As in previous sections, we explore (177|) in the two limiting regions: far from (r — )■ +oo) and 
close to (r — t- 0) the brane. Taking into account (121 p . in both regions the equation reduces to 

- + 9J«o(r) = 0. (78) 



2 

So, the solution in these to regions has forms 

«o(r)L,^o = Cle-W^ «o(r)|,^+^ = C^e'^^/^ , (79) 
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with Ci and C2 being some real constants. Then, using f l2T]) and fl79|) . in these regions the zero 



mode time-averaged Lagrangian Lfj^ 



-g (^'^iT^DA'^) has the following forms 



2 



r=0 2 



(0) 



RoC2e-'''iPix'')iYduHx>' 



^0) 



On the brane the 6D Lagrangian has the standard 4D form, and the integral over extra coor- 
dinates r and 6 in in the action ( l69l) is convergent. So the 4D fermions are localized on the 
brane. 



6 Localization of Gravitons 

For spin 2 gravitons we consider the following metric fluctuations 

ds^ = e'^"' {g^^ + V) dx^'dx'' - dr^ - Ry^'-^''de^ , (81) 

where 

g^, = (1, -e", -e", -e") . (82) 
It is straightforward to show that equations of motion of the fluctuations h^i, are 

{V^ g^'^'dNKJ) = . (83) 

These equations of motion for the fluctuations in the present background are equivalent to 
that of scalar field f l24p considered in Sec. 3. Accordingly, the localization problems for spin 2 
graviton and scalar field are similar, and gravitons are also localized on the brane. 



7 Summary and Conclusions 

In the article we have introduced new nonstationary 6D standing wave braneworld model 
generated by gravity coupled to a phantom-like scalar field and have explicitly shown that 
spin 0, 1, 1/2 and 2 fields are localized on the brane by the universal and purely gravitational 
trapping mechanism. In our model, as opposed to earlier static approaches with decreasing 
warp factors [imfTT]. localization takes place in the case of increasing warp factor. 

Here we have considered the model for u/a = Zi, where Zi is the first zero of the Bessel 
function Ys. But according to condition ( ITSl) there is possibility to control the number of "zero 
circle" islands in the 2D extra space that can provide interesting approach to the old problem 
- the nature of flavor. 
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